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Abstract 

Spin motive force induced by the Rashba interaction in the presence of strong sd interaction be- 
tween conduction electron and localized spin is theoretically studied. The motive force is calculated 
by evaluating the time-derivative of the current density on the basis of microscopic formalism. It is 
shown that there are two motive forces, one proportional to x n, the other, perpendicular com- 
ponent, proportional to x (n x n), where and n are the Rashba electric field and localized 
spin direction, respectively. The second type arises in the strong sd coupling regime from the spin 
relaxation. The appearance of perpendicular component from the spin relaxation is understood 
from the analogy with the current-induced torques. In the case of domain wall motion, the two 
contributions to the spin motive force are the same order of magnitude, while the first term dom- 
inates in the case of precession of uniform magnetization. Our result explains appearance of the 
perpendicular component in the weak sd coupling limit, recently discussed in the context of spin 
damping monopole. Detection of AC voltage induced by the precession of uniform magnetization 
serves as a experimental evidence of the Rashba interaction in films and wires. 
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I. INTRODUCTION 



In electronics, electric motive force driving electron's charge is used. Recently, control 
of electron's spin has been carried out in spintronics technology, and use of another motive 
force, spin motive force, has become possible. The idea of spin motive force induced by 
an sd interaction between conduction electron and localized spin was argued by Berger in 
the case of a domain wall motion in 1986 In the strong sd coupling regime, the spin 
motive force is described by a U(l) gauge theory for an adiabatic component, Al „, of the 
spin gauge field, (// = x, y, z, t and a = x, y, z are indices for space time and spin) ^|. 
In the adiabatic limit and in the absence of spin relaxation, Al has a U(l) gauge symmetry, 
and the spin electric field, defined by spin motive force divided by the electron charge (e), 
is E Ki = dtAli — diAl p[ [J. By use of a unit vector, n(r,t), describing the local spin 



direction, it is E B) i = —7^n ■ {n x c^n)^ This spin electric field was experimentally observed 



in the case of moving domain wall 

When there is spin relaxation, the spin motive force is modified. In the context of 
current-driven torque, spin relaxation induces a perpendicular component, /3 sr (n x (j ■ V)n), 
to the adiabatic spin-transfer torque (j ■ V)n (j is the applied current density and /3 sr is 
a coefficient proportional to the spin relaxation rate) jf], [7\. As discussed by Duine jg] 
and Shibata (9], spin motive force is the inverse effect of current-induced torque. The 



spin relaxation is thus expected to modi 
E^ oc (3 sr n • [n x (n x V,n)] = h ■ Vjn 



spin electric field to give rise to a correction, 
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The above results apply irrespective of the origin of the spin relaxation; Relaxation due 
to spin flip scattering by magnetic defects and that due to spin-orbit interaction at heavy 



impurity sites result in the essentially the same result 12|,ll3|. In contrast, Rashba spin-orbit 



interaction, which arises from a breaking of the inversion symmetry, leads to a qualitatively 



and to give rise to a spin Hall effect 
electric field was studied very recently 
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-|l6[. The effect of Rashba interaction on the spin 



regime and found a spin electric field of i? s = a[JE_oX (n x n)] = E s ±, where a is a constant 
and .Er is the electric field of Rashba interaction 171 ] . They also calculated the spin magnetic 
field, B s , and found that the spin electromagnetic field has a monopole. Namely, the field 
has a finite monopole density and monopole current density, defined in electromagnetism 



theory as p m = V • B s and j m = V X E 8 + B s , respectively. The method employed by 
Takeuchi to address the spin electromagnetic field is unique; the fields were identified by 
calculating the induced electric current density and then comparing the result with a general 
expression, j = aE s + V x B s — DVn e , where a is conductivity, D is diffusion constant 
and n e is charge density. This approach is highly useful to study the weak coupling regime, 
where adiabatic component of spin gauge field can not be defined. Kim et al. studied, on 
the other hand, a strong sd coupling regime and obtained a different form of E s = b[En x h] 
(b is a constmt , Q. Therefore, l e spin electric fields so far identified are different for 
strong and weak coupling regimes. The absence of E s ± in the strong coupling limit would 
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be due to the fact that the effect of spin relaxation is not taken account of in Ref 
is suggested from the case of current-induced torques. 

The aim of this paper is to study the spin electric field (motive force) induced by Rashba 
interaction and sd interaction on a basis of a microscopic theory taking account of the spin 
relaxation due to spin-orbit interaction driven by heavy impurities. Spin relaxation effect 
is expected to be essential in grasping global behavior of the spin electric field. We have 
mentioned above two methods used for far to calculate effective spin electromagnetic field; 
1) by use of adiabatic gauge field, and 2) by calculating electric current. Here we will use 
different approach. The motive force, F, acting on conduction electrons is defined by an 
average of the time derivative of electron velocity multiplied by mass, m. It is therefore 
calculated by evaluating ^ as 

m dj 



en e dt 



(1) 



where n e is the electron density and e is the electron charge. This argument was employed 



in Ref. 18] for a phenomenological derivation of spin motive force. We will estimate Eq. 
([1]) microscopically by rewriting the time-derivative by a commutator of the current density 
operator, j, and the Hamiltonian H as 

' m '[H,J]), (2) 



ehn 



where ( ) denotes quantum average. We will show that when spin relaxation is taken into 
account, E s> ± indeed emerges even in the strong, sd coupling regime, as is consistent with 
previous analyses on current- induced torques jslQ, [l2, 13] and on spin motive force without 



Rashba interaction 



a. 



In fact, the Rashba interaction gives rise to a spin electric field 



E S R = 7r [(a R x n) + /3r(cxr x (n x n))] , 



(3) 



where 7r = ^ and represents the strength of the spin relaxation. The first term was 
pointed out in Ref. 



18j, and the second term has the form found in Ref. in the weak 
sd coupling limit. In the case of precession of uniform magnetization, the first contribution 
dominates, since /3r is small (of the order of 0.01), as discussed in Sec. IIVI For a moving do- 
main wall in a perpendicular anisotropy system, the contributions from the first and second 
contributions of Eq. are the same order. For current-induced magnetization switching 
in a perpendicular magnetic random access memory, the first contribution enhances the 
damping and consequently, increases the threshold switching current density. 

Based on a series of experiments on Pt/Co/A10 x systems 20|,|2l|, Miron et al. argued that 
an injection of an in-plane current into this system generates a large transverse magnetic 
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field, predicted from Rashba theories Q, [15|. On the other hand, Liu et al. 
reported that spin transfer torques induced by an in-plane current in similar systems can 
be quantitatively explained by the spin Hall effect and they could not find any signatures of 
Rashba effect in their measurements. These contradictory experiments have led to extensive 
discussions on the existence of Rashba effect in ferromagnets. 

As we will show in Sec. IIV[ the amplitude of Rashba-induced alternating spin voltage 
in the case of a precession of a single domain magnetization of a size of lOOnm is predicted 
to be as large as 0.2mV at a frequency of 100MHz if the Rashba coupling constant is 
3 x 10- 10 eV-m 



«R 



. From the Rashba spin motive force, Eq. fl3]), the Rashba coupling 
constant, «r, is determined since the prefactor, 7r, is universal and does not depends on 
material properties is small and is neglected in a uniform magnetization case). At the 



same time, ckr is measured from a domain wall speed |20|. Measurement of AC spin voltage 
without domain walls is thus expected to provide a highly useful information concerning the 
existence of Rashba effect. 

For confirmation of Rashba-induced spin motive force, it is essential to exclude similar 



signa! from the spin pumping and inverse spin Hall effects & In the case of perpendicular 
easy axis systems like in Ref. [20[, the dominating part of inverse spin Hall signal vanishes, 
since the polarization of spin current is parallel to the spin current flow. 

In the next two sections, we show derivation of the Rashba-induced spin motive force. 



Those who are interested in the result and its physical consequences only may skip the 
section and read Sec. [IVl 



II. SPIN MOTIVE FORCE WITHOUT SPIN RELAXATION 



18] by evaluating Eq. 



We first study the case without spin relaxation, studied in Ref. 
([2]). The interactions we consider is the sd interaction and the Rashba interaction. The 
Hamiltonian is 

h 2 v 2 



H 



d rc 



2m 



e F - A sd n ■ a + ict R ■ (V x <x) 



(4) 



where A s d is the energy splitting of the conduction electron, n represents the direction of 
localized spin, ckr is a vector representing direction and strength of the Rashba interaction. 
The electron is described by creation and annihilation operators, c' and c, respectively. 
These operators have two spin components as c = (c + ,c_) where ± represents the spin 
index. The electric current in the present system reads 

e 



c T V c 

2m 



h 



cUolr x er)c. 



(5) 



We consider a strong sd coupling limit, and carry out a local gauge transformation defining 
new electron operator by a = Uc, where U is a 2 x 2 matrix. The matrix is chosen to 
diagonalized the sd interaction, i.e., U = m ■ a, where m = (sin | cos0, sin | sin0, cos |) (9 
and (f) are the polar coordinates of n) 25j. After the gauge transformation, the Hamiltonian 
and the current read 

h 2 v 2 

- e F - A sd a a 




2m 



2 Al 



+hA[ a [ aia + - 



3 + e ijk a Ki R k i ) a f 

1 fh 2 (A s ) 2 



m 



— tiikOtRiRkiA, 



a} a 



(6) 



where At 



-J \ti[aiU^ dfJJ] = (m x d^m)i is the spin gauge field and Rm = 2m k mi — 5m is 
rotation matrix element. 



By calculating the commutator [H,j] in Eq. (j2]), we obtain the force operator as 
- _^ c t V (n ■ <r)c + ^S^ct[a R x (n x <t)]c + 0((a R ) 2 ). 



h 2 



(7) 
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FIG. 1. The Feynman diagrams describing the spin motive force force arising from the Rashba 
interaction at the lowest order. The vertex with E R represents the force induced by the Rashba 
field (the second term of Eq. fl7D), and the interaction with the spin gauge field, A s , is included to 
the linear order. 



The first term containing Vn is the term discussed in Refs. 26|, |27| in the context of 
current- driven domain wall motion. Carrying out the gauge transformation and taking the 
field expectation value, the force reads 



Fi = — Y] I (Win^Rjn - ^ V] e ijk e k i m a Rj niR mn \ (aV n a) + 0(a R 2 ). (8) 



n. 



Y I (Vinj)Rjn ~^Y e 

jn \ klm 



The contribution from the first term of Eq. (jHJ) turns out to be the second order in A B , 
and is neglected in the adiabatic limit. The contribution from the second term, described 
in Fig. [U is calculated as 

F i = ~i ~ S Y Yl e ij^klmOiRjniRmn Y 4" Y e~ i<? V™ Y tT K#fc^a£fc +( ^ + d < , 
6 jn klm a qfi koj 

(9) 

where g kul is the electron Green's function on the Keldysh contour with wave vector k and 
angular frequency u, < denotes lesser component and = J We have neglected the 
contribution containing the spatial component of the spin gauge field, since they turn out 
to be order of (AJ 2 . The lesser component is calculated as 

(10) 

Trace in the spin space is taken by use of 

tr[a n Aa a B] = (S na - S nz 5 ao ) Y a ± b t + $n Z 5 ao Y A ± B ± ~ ie naz ^(±)A ±J B T , (11) 

± ± ± 
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where A 

given by 

_ 2mA sd 
hn e 

where 



401 




' B + \ 


and B = | 


I A 1 




l 5 J 



are diagonal matrices. The force is therefore 



jn fcZm \ a / 



(12) 



-q,W+fi, J 

+n,±] 



7i — ih'y] [ fu+n9koj±9k+ q ,u+n,T + /wg/L±fi , fc+q,w+n,=F + (fu+n fu)9kw±9k- 

12 = [ — ^+ n 5'fcw±5 , fc+q,w+^,± + fu9ku>±9k+q,u)+n,± + (/w+n ~~ fu)9kuj±9k- 

kui± 

13 — ~^'^^(^ z ) [~ f^+^9ku)±9k+q,ui+n,T + f^9kui±9k+q,ui+n,T + (/w+fi — /w)#Lj±.S , fc+<j,w+n,=F] • 

(13) 



The terms proportional to 7° and 7° are calculated using 25 ] 

^ ^ ^klm^ljRmnAg Q -R mz v4 s g) = 
Iran 

and 

1 <9n. 



2 dt 



(14) 



2 at 



(15) 



and the term proportional to 7° vanishes. The result of the force is thus 



tin 



[7i(atR x n) + 7°(a R x (n x n))] 



(16) 



We evaluate the coefficient 7° and 7° in the clean limit (infinite elastic lifetime) and in 
the adiabatic limit (approximating q = in the Green's functions) using the identity 



The result is 



7?- 



2A 



9L± - 9L± = 27iiS{hw - e k ±). 



(17) 



k± 



■(n+ -ra_), 



(18) 



where n± = J2k f( e k±) is the density of spin ± electron, and 

^ E (/(*«=) - /(*f )) = °" ( 19 ) 



7 3 °- 



fe± 



The motive force on the electron spin in the absence of spin relaxation is therefore ob- 
tained as 

F = ~(a a xn), (20) 
where s = n + — n_ is the spin polarization of the conduction electron. (In the strong 



coupling limit, s = n c .) The perpendicular component (ckr x (n x n)^ 
in the absence of spin relaxation, in agreement with the result in Ref. 



thus does not exist 



III. EFFECT OF SPIN RELAXATION 



Here we study the effect of spin relaxation on the spin motive force. The spin relaxation 
we consider is the one due to spin-orbit interaction due to heavy impurities, described by 
the following Hamiltonian: 



J d 3 r^2e ijk (V l v£ ) )c j Vj a k c, 

ijk 



(21) 



where vio is random potential due to random impurities. To obtain physical force, the 
potential is averaged over the random distribution, and the averaging is carried out taking 
account of the spin direction, denoted by the label k as carried out in Ref. 13]. After the 
gauge transformation, the spin-orbit interaction reads 



/ d 3 r eijk(Viv!$)Rki ( -^« f Vj a x a + A l sj a ] a 

ijk 



(22) 



and, in the Fourier representation, 

H so = -i E E e v* v so\ k ' ~ k - q)(k' - k - q)i ^ \( k> + k)jR k i{q, Q)a\ u aia k/> , 



kk'q ijk 



q'Q,' 



(23) 



Here Fourier components are defined as A\Aq,Q) = f d 3 r f dte tq r e %nt A l S j and v so (q) = 
J d 3 re iq " r v B0 (r). We assume that the impurities are heavy and thus the momentum k' and 
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FIG. 2. The Feynman diagrams describing the corrections of the Rashba- induced spin motive 
force due to the spin relaxation, F^ ,A . These contributions turn out to be higher order in the 
spin gauge field, and thus neglected compared with the contribution of Fig. [3j 

k in Eq. (1231) can be treated as independent. In contrast, the momentum and the frequency 
carried by the spin gauge field and rotation matrix, q, q', Q and Q', are small compared to 
k and k' considering an adiabatic regime. The Hamiltonian of spin-orbit interaction then 
reads 



on 



(24) 



kk'q ijk 

q'fl' 

Let us proceed to estimate the correction of the spin motive force due to the random 
spin-orbit interaction. Since the bare spin motive force without relaxation is described by 
the Feynman diagram shown in Fig. one would guess that the relaxation correction is 
represented by the processes shown in Fig. [21 These processes, however turn out to vanish 
at the first order in the gauge field as we will now demonstrate. 

The sum of the three contributions in Fig. [2], denned as F^' A , in the adiabatic limit 
(q ->■ 0) is 

r-,(sr).A 2mA sd 

F i = — 

hn c 

ijklmn afi-/ 

X ^defhkfk'XRaaWRdpW (vg> (k! ~ k)vg> (-k' + k)) 

abcdef 

x X^ tr [ (Xn 9kuj ''y9k,uj+tt a 'a9k>,uj+n 1 +n a a 'i39k,uj+n 

+ < J n9kuj 'a9k>,Lu+n 1 cr i39k,Lu+n 1 +n2 a 'y9k,w+n 
+ a ri9ku) a a9k>,w+n 1 a 79k',u)+n 1 +n 3 cr l39k,u>+ri\ > 
where f2 = + f2 2 + ^3- The average of the random potential is calculated as [1 

(vi:\k>-k)vif(-k> + k)) = (vi)5 ad 



f2i ,£^2,^3 



kk' 



(25) 



(26) 
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where (v 2 Q ) represents the squared average of the spin-orbit strength. Using the rotational 
symmetry, we replace k^kf — > ^k 2 Sbf in Eq. f l25|) . The trace in the spin space is calculated 
using the following identities. 

° a g(T a = 9 + ^9 

a 

a c t 9cr y fcr a = 25 yz a z (gf - gj) - ga^f, (27) 

a 

where g = go + o z 9\ is an y diagonal matrix, ~g = g — o z g\ and similarly for /. The result is 

( 7 { sr) ' A (a R x n) + 7? r) ' A («R x (n x n))) , (28) 



,(sr),A _ "^A srf / ( sr ),A, 



1 hn e 
where 

7 ( sr )' A = (v 2 so ) I p- ^{±)k 2 {k'f [f u +n\gL±9k^n&(9u fi ,+n,T + 2 9k'< u +si,±)9k< u +ti& 

,uj,±9k,u+Q,^f 9kuj±9k' ,u,±9k' ,Lu+Q,^9k,uj+Q,^f\ full? ^ a J 

+ (f u +n ~ fJ)\gL±9k ,uj+n,^f(9k',uj+n,^f 2 9k',u)+n,±)9k,uj+n,T 

+9k (9 ] , (29) 

and 7| sr ^'" 4 is defined without the sign (±) in the summation. The term 'y^^ gives only 
a small correction to Eq. ( 120]) . and we are here interested only in , y^' A . In the adiabatic 
limit Q — > 0, the leading contribution turns out to be 



73 ^ — 2- / — (±) (An-n /w)[25'fc a; -|-(5 , fc itJi = F ) 2 5 , fc' iWi -|- #L;T(#fc,w,±) 2 #fc',w,±] 
71 fcfc'± 

= O(fi). (30) 



The contribution of 73 Sr ^ ,j4 to a force is therefore a second order in the time derivative as 
a n x J^(n x n). It is neglected in the adiabatic regime, and the processes shown in Fig. [2] 
do not thus to contribute to a perpendicular spin electric field proportional to ckr, x (n x n). 

We will now show that dominant contribution arises from simpler processes shown in Fig. 
|3l In fact, in those processes, the rotation matrices, Rij, at the two vertices of the random 
spin-orbit interaction are defined in general at different space and time coordinates. At equal 
time, Rijif)Rjkit) = 5ik, and this relation results in a vanishing contribution to the force, 
while non-trivial contribution is expected at different time, since ^ Rij(t)Rjk(t') ^ 5ik if t ^ 

10 
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FIG. 3. The Feynman diagrams describing the spin motive force arising from the Rashba inter- 
action (En) and spin relaxation, represented by the random spin-orbit interaction, v so . Although 
there is no explicit interaction with spin gauge field (A s ) in the diagram, the contribution turns 



out to results in a force containing the spin gauge field, proportional to qr x (n x n) 



28|. 



t'. (In the frequency notation, R^ carries a finite frequency, and thus ^ Rij(fli)Rjk(fi2) ^ 
5ik if fii 7^ ^2-) This fact in fact leads to a perpendicular contribution as we will now 
demonstrate. The contribution of Fig. [3J F^ ST \ after averaging over random potential 
according to Eq. (1261) is 



( sr ) AmA sd \ ^ \ ^ \^ intsr^ 

1 ~ ~ % ~§~hn~ ^ 2^ €i i k€klmaR i niKm -t 1^ e 1^ 

6 ijklmn afl-y f2i,f!2 fcfc' 

Xk^k'fRn^Rn^)^) 

] , (31) 

where f2 = + f2 2 - It is easy to see that the contribution from the trace evaluated at 
S7i = SI2 = vanishes, and the leading contribution arises from those at the linear order. 
Summing the contribution obtained by replasing a with and Qi with Q2, we can rewrite 
Eq. (ED as 

F i ST) = "^g - ^^ £ S e ijk e k i m a Kj riiR mi ^ e*™ ^ 

6 ijklmn a/37 f2i,f22 fefe' 

x J^tr [(7 7 g ku _a(7 a g k , tUJ+An a^g kiU}+ n +M fcu -§ (7 ^fc>-^3ik, u+ n_ > ( 32 ) 

where Af2 = — ^2)- We note here that the contribution to Eq. (132)) linear in f2(= 
Qi + Q 2 ) vanishes because of ^2 n R na (t)R n i3(t) = S a p. We can thus expand Eq. (132|) with 
respect to Afl at f2 = 0. We also know that the contribution asymmetric with respect to a 

11 



and (3 contributes to the force. The result of the force is therefore 



F 



(sr) 



.4 mA, 



-%— 



9 hn P 



^2 ^ e vk^ki m aiijniR mi ^ e jni (fii-fi 2 )^ 



ijklmn a/3"/ 



ni,n 2 



kk' 



X 



k^k'fRU^R^) (v 2 so ) h 



doo 
2^ 



x [f'(u)ti [cr^ciaigl,^ - gi,Jopg k ^ - (a /?)] 
The summations over Qi and f2 2 are carried out as 



(33) 



/] e in \Qi — Q 2 )Rna(^l)Rnl3(^2) = Rna(t) Rnfiit) — R na (t) R n/3 (t) . (34) 



Carrying out the summation over the index n, 



y^;{R na (t)R n p{t) ~ R na (t)Rnp{t)) = ^{m a m p - m a rr 



(35) 



Trace in Eq. ( 133]) gives rise to terms proportional to e a p z 5 lz , tap-y and to 5 ai 5p z — 5^5 C 
The contribution containing 5 1Z vanishes, since it is proportional to e kim n iRmz 
J2im e kim,nin m = 0. The contribution containing e Q( g 7 leads to (by use of Eq. ( |34|) ) 



^niR mz {t) R na (Qi) R n/3 (fl 2 ) 
A dn k 



I III 



[mf.m z + m k m z 



dt 



Similarly, the contribution containing 8 ai 5p z — 5p 7 6 az leads to 

^klrn 

(Sa-ydpz — S/3y5 az )niR mz [t)R na {Qi)R n p(Q,2) 
- 4l " x f) t - 



Defining 



16 



/dcu 
— tr [/'(a;) [<r 7 <4,<r a (<& )£ , - sj^Jtr^ - (a o /?)] 

y (^) 



>r) 



(sr). 



(36) 



(37) 



(38) 
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the force is (the contribution from the term vanishes) 

F * (sr) = iif (7!sr)(c * R x h) + 7 * sr)(aR x (n x (39) 



Here again we focus on the term 73 , which reads 

7? r) = ~ (v 2 so ) k 2 (ky(g^ ± - gl,, ± W k± gl, T - <?^, ± ). (40) 

Summation over are carried out as 

5>T(<&, ± - 9l>,±) = 2tt^ ± (A; f± ) 2 
fe' 

E fc X±^,T = ^M** + ) 2 + M^f-) 2 ), (41) 

where and fcp± are the density of states and Fermi wave vector of the conduction electron 
with spin ±, respectively. The result is thus 

7? r) = «o> (M** + ) 2 + »-(k F -) 2 ) 2 . (42) 

The same diagram as in Fig. [3] was discussed in the context of the Gilbert damping in Ref . 

(sr) 

in the case of relaxation due to random magnetic impurities. The result of 73 is in 



agreement with their result. 

IV. APPLICATION 

In this section, we apply our results of spin motive force. We investigate three cases, a 
precessing uniform magnetization, domain wall motion, and current-induced magnetization 
switching, to explore the effect of spin electric field in detail. Adding the contribution from 
the adiabatic gauge field {2, 8], the total spin electric field is written as 

E S j = 7 [n • (n x V l n) + /3 7 n ■ V»n] + 7r [(a R x n) + /3 R (a R x (n x h))} . , (43) 

where 7 and 7r represent coupling to the adiabatic gauge field and to the Rashba field, 
respectively. In the strong coupling limit, 7 = ^ and 7r = ^- Coefficients /3 7 and /3r 
represent the strength of the spin relaxation on spin electric field. These constants are of 
the same order as the /3 coefficient of the spin-transfer torque. 
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A. Precession of uniform magnetization 



We first consider the spin electric field driven by a precessin g u niform magnetization. 



The effect of perpendicular component, E s± , was studied in Ref. [17J. The adiabatic gauge 
field contribution, the term proportional to 7 in Eq. ( I43p . vanishes identically for a uniform 
magnetization. We consider a film or a thin wire with a perpendicular easy-axis anisotropy 
along z axis. The Rashba field is also along z axis (ckr = or!) 20(. We consider the 
magnetization precessing with a constant tilt angle 9. The direction is represented by a unit 
vector 



n = (sin 9 cos 4>(t), sin#sin</>(i), cos#). 
Its time derivative is n = sin^e^, and the result of spin electric field, E s = E, 



(44) 



R 



^ COS(0 + if) ' 



-7r«r0 sin 9 + (/3 R ) 2 cos 2 9 



sm(cj) + if) 




(45) 



where f = sin 1 



/3r, cos 9 



--. In the uniform magnetization case, contribution from spin 



v/1+OSr) 2 cos 2 ' 

relaxation (-E s ± proportional to /3r(^C 1)) is small compared with E s «. The magnitude of the 
electric field is thus estimated to be \Ef-\ = 7rO!ru; sin 9, where u is the angular frequency of 
precession. Considering a case of strong Rashba interaction, «n = 3 x 10~ 10 eV-m 18[, the 
field is estimated for a frequency of v — u/(2ir) = 100MHz as \Ef-\ = 2 x 10 3 V/m. For a film 
with size of lOOnm, the voltage expected across the film is 0.2mV. By detecting this AC spin 
motive force, therefore, one can verify the existence of the Rashba interaction and estimate 
its strength in films. A great advantage of using spin motive force is that the voltage signal 
accumulates over the system as far as it is single domain, in contrast to the domain wall 
case, where the signal arises only from inside the domain wall. Rashba-induced spin motive 
force would thus be very useful to distinguish the Rashba effect on the spin-transfer torque 



22 



23 



29) 



and spin Hall effect, which is an urgent issue at present 

A DC component of voltage signal generated by a magnetization precession in a sys- 
tem of Pt film in contact with a permalloy film was observed in Ref. 24[, and the result 
was explained based on a spin pumping and inverse spin Hall effects. The AC component 
of the voltage may contain Rashba-induced motive force if there is Rashba interaction at 
the permalloy-Pt interface. Bilayer or multilayer systems with a perpendicular easy axis 
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anisotropy like in Ref. 20( is suitable for detection of Rashba-induced spin motive force, 
since the effect of inverse spin Hall effect is neglected. This is because the spin polarization 
a of spin current pumped by the magnetization is parallel to the spin current flow j s and 



thus inverse spin signal, proportional to j s x a 



24l | . vanishes. 



B. Moving domain wall 

Let us consider the effect of the spin electric field on moving Bloch domain wall in a thin 
film with a perpendicular easy axis anisotropy. The film is in the xy plane, and so the easy 
axis is along z direction. The Rashba field, = E-&Z, is chosen along z direction, as would 



be the case in experiments 
coordinates as 



20j. The magnetization unit vector is written by use of polar 



n = (sin 9 cos 0, sin 9 sin 0, cos 9) . 



(46) 



A planar domain wall with magnetization changing along x direction and at position x 



X{t) is described by a solution sin 9 



i 



cosh ; 



25] . where A is wall thickness, treating as 



- s -^e e and 



uniform inside the wall, we thus obtain V x n ■ 

h = sin 9 ^0e + y e^j , (47) 

where = (— sin 0, cos 0, 0) and eg = (cos 9 cos </>, cos 9 sin <p, — sin 9). The total spin electric 
field averaged inside the domain wall (using (sin^) = - and (cos 9) = 0) reads 

M 





v / 



^ cos (b ' 



U+R X 



sm< 




(48) 



/ 



The first contribution without Rashba interaction was discussed previously by Lucassen et 



al. 



We see that the contribution proportional to wall speed, X, arises solely from 



spin relaxation (contributions containing f3 or /3-g), while contribution arises without spin 



271 ]. the Rashba 



relaxation. In the case of a steady spin-torque motion with oscillating 
contribution, proportional to cos or sin leads to AC component as predicted by Kim et 
al. isi, while the adiabatic gauge field contribution (7 contribution) has a DC component 
if and X have finite average. Experimental detection of Rashba-driven AC contribution 
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rotating in the xy-pl&ne provides a novel mean to confirm the dynamics of the wall's internal 
degrees of freedom, 0, after the Walker's breakdown in the field-driven case and in the 
current-driven motion in the intrinsic pinning regime 27], |30 1 . 

Let us study our result, Eq. (jUJ), in more detail considering the case without extrinsic 



pinning potential. As described in Ref. 25|, the equation of motion in this case is (neglecting 
oder of a 2 , a/3 <C 1, where a is the Gilbert damping constant) 

dX 

(49) 



^ = ^{{(3-Pa)j-asm2 

where P is spin polarization of the current, v c = + (K± is hard axis anisotropy energy 
and S is the magnitude of spin). 

The Rashba contribution to the spin electric field is governed by a factor of (neglecting 

0(P 2 )) 



+ Aiyj = y {{13 + P(/3r - a))j + (fo - a) sin 2 
The Rashba contribution to the spin electric field is therefore 



(50) 



^ cos 6 » 



E, 



R = -^E R ^- [((3 + P((3 R - a)) j + 2((3 R - a) sin cos 0] 

7T A 



smi 




(51) 



/ 



By observing the AC component, parameters (/3r — a) and (/3 + P{(3 R — a)) are accessible 
experimentally. Note that the spin relation contribution proportional to /3 R (the fourth term 
of Eq. (|43p ) contributes to the same order of magnitude as the bare contribution (the third 
term of Eq. (|43l) ). since X / A is larger than by a factor of 25]. 



C. Threshold current for magnetization switching of a perpendicular layer 

For a single domain state, the total spin electric field (Eq. (115]) ) is reduced to 



Ef = ±^[{z x h) + (3 R (z x (n x n))]. 



(52) 



This spin electric field generates a spin current density via j s = cr^E+ — cr^E s = cr c E+ 
where ovn is the electrical conductivity of majority (minority) electrons and a c = <r^ + a±. 
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The spin transfer torque due to? s can be o 
and nonadiabatic spin torques [14j, |l5j, |29|, [31 



atained by inserting j s into Rashba adiabatic 
, |32| (neglecting order of /3r/3) 



T{js) = lg°c (^) [nx{D-h) + (3 R n x (f x (D • n)) - (3n x n x (D ■ h)\, (53) 

where 7 9 is the gyromagnetic ratio, and D is a 3 by 3 tensor with the element Dy = 
Sij(l — 5i Z ). Then magnetization dynamics with no current injection is described by the 
modified Landau-Lifshitz-Gilbert equation, 

h = —~f g n x H e ff + an x h + an x (D ■ n) +a[/3 R n x (fx (D-h)) — (3n x n x (D -n)], (54) 

where i? e // is the effective magnetic field acting on n and a = 'y g a c (maR/eh)' 2 . One finds 
that the third term on R.H.S. of Eq. (15"4"1) has the same form with the Gilbert damping 

n 

torque and thus makes the damping anisotropic [1 81 ] . 

We next discuss the role of the last term on the R.H.S. of Eq. ( 15~4"|) . which is related to (3 R 
and p. For a perpendicularly magnetized layer where magnetization undergoes a small-angle 
precession around the easy axis, n and -tf e // are described by (m cos (cut), mo sin (cut), 1) and 
(—N xy M s mo cos(u)t), —N X yM s rriQ sin(a;£), Hk — N Z M S ), respectivley, where mo C 1, w = 271"/, 
/ is the precession frequency, and N xy (N z ) is the in-plane (out-of-plane) demagnetization 
factor. For this case, the precession frequency / is given by / = ^ Hk jj^^z^ff^ and 
thus the last term on the R.H.S. of Eq. (|54|) may be interpreted as a modification of the 
gyromagnetic ratio induced by spin electric field. 

Finally we investigate effect of Rashba spin-orbit coupling on current-induced magneti- 
zation switching in perpendicular spin valves, which is of considerable interest for ultrahigh 
density spin transfer torque magnetic random access memories (STT-MRAMs). To find a 
threshold current for perpendicular STT-MRAM of which free layer is subject to Rashba 
spin-orbit coupling, one uses Eq. f )54|) with including the Slonczewski's spin torque term, 
7 g 2eMgt F n x ( n x z), where rj is the spin torque efficiency factor, j e is the charge current 
density, Ms is the saturation magnetization, and ip is the thickness of free layer. Assuming 
a small-angle precession (i.e., mo 1) and using the energy dissipation rate (dE/dt) (i.e., 
(dE/dt) = —M s (H e ff ' dn/dt) = 0, where (...) is a time-average over a period), one finds 
a threshold current density j S w as 

17 



Jsw = ( a + a)-^^-L(H k - (N z - N xy )M s ). (55) 

From Eq. (|55[) . one finds that Rashba spin-orbit coupling enhances the damping by a and 
consequently increases jsw- On the other hand, the modified gyromagnetic ratio does not 
affect jsw because jsw is entirely determined by the competition between the damping 
torque and the Slonczewski's spin torque. One may want to reduce jsw as much as possible 
for application of perpendicular STT-MRAM. Eq. ( )55l) suggests that Rashba spin-orbit 
coupling should be minimized for this purpose. 



V. SUMMARY 

We have calculated the spin motive force acting on the electron induced by the Rashba 
interaction and sd interaction based on a microscopic formulation. We have shown that there 
are contributions in the limit of strong sd interaction; one proportional to .Er x h(= E s ^), 
which arises without spin relation, and the other perpendicular component proportional to 
x (n x n)(= E s ^) arising from the spin relaxation (E^ and n are the Rashba electric 
field and localized spin direction, respectively). 

Rashba-induced field arises even in the uniform magnetization case, in contrast to the 
spin electric field arising from adiabatic gauge field, which is finite only if magnetization has 
both spatial structure (finite Vn) and time-dependence. Measurement of spin motive force 
driven by uniform magnetization precession would be useful in confirming the existence of 
the Rashba interaction in multilayer systems with perpendicular easy axis like the system 



studied by Miron [20] 



In view of our result, the appearance of a perpendicular spin motive force in the weak sd 



coupling limit discussed in the context of spin damping monopole [17| is naturally under- 
stood, since spin relaxation is strong in the weak coupling regime having no particular spin 
quantization axis. There is a possibility that the other spin electric field E<q also contains 
monopoles. This possibility is currently under study by investigating the spin magnetic field 
in the strong coupling limit. 
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